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This Letter studies the decoherence in a system of two antiferromagnetically coupled spins that interact 
with a spin bath environment. Systems are considered that range from the rotationally invariant to highly 
anisotropic spin models, have different topologies and vàlues of parameters that are fixed or are allowed to 
fluctuate randomly. We explore the conditions under which the two-spin system clearly shows an evolution 
from the initial spin-up - spin-down state towards the maximally entangled singlet state. We demonstrate 
that frustration and, especially, glassiness of the spin environment strongly enhances the decoherence of the 
two-spin system. 
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The interaction between a quantum system, called 
central system in what follows, and its environment af- 
fects the state of the former. Intuitively, we expect that 
by turning on the interaction with the environment, the 
fluctuations in the environment will lead to a reduction 
of the coherence in the central system. This process is 
called decoherence G] ■ In general, there are two dif- 
ferent mechanisms that contribute to decoherence. If 
the environment is dissipative (or coupled to a dissi- 
pative system), the total energy is not conserved and 
the central system + environment relax to a stationary 
cquilibrium state, for instance the thermal equilibrium 
state. In this paper, we exclude this class of dissipa- 
tive processes and restrict ourselves to closed quantum 
systems in which a small, central system is brought in 
contact with a larger quantum system that is prepared 
in its ground state. Then, decoherence is solely due to 
fact that the initial product state (wave function of the 
central system times wave function of the environment) 
evolves into an entangled state of the whole system. The 
interaction with the environment causes the initial purc 
state of the central system to evolve into a mixed state, 
described by a reduced density matrix[3], obtained by 
tracing out all the degrees of freedom of the environ- 
ment pmiiE]. 

Not all initial states are equally sensitive to decoher- 
ence. The class of states that is "robust" with respect to 
the interaction with the environment are called pointer 
states %. If the Hamiltonian of the central system is 
a perturbation, relative to the interaction Hamiltonian 
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H int , the pointer states are eigenstates of H int E] . 
In this case, the pointer states are essentially "classi- 
cal states", such as states with definite particle posi- 
tions or with definite spin directions of individual par- 
tides for magnètic systems. In general, these states, 
being a product of states of individual partides form- 
ing the system, are not entangled. On the other hand, 
decoherence does not necessarily imply that the cen- 
tral system evolves to a classical-like state. If Hi nt 
is much smaller than the typical energy differences in 
the central system, the pointer states are eigenstates of 
the latter, that is, they may be "quantum" states such 
as standing waves, stationary electron states in àtoms, 
tunnclling-split states for a particle distributed between 
several potential wells, singlet or triplet states for mag- 
nètic systems, etc. [ü]. This may explain, for example, 
that one can observe linear atòmic spectra - the initial 
states of an atom under the equilibrium conditions are 
eigenstates of its Hamiltonian and not arbitrary super- 
positions thereof. 

Let us now consider a central system for which the 
ground state is a maximally entangled state, such as a 
singlet. In the absence of dissipation and for an envi- 
ronment that is in the ground state before we bring it 
in contact with this central system, the loss of phase co- 
herence induces one of following qualitatively different 
types of behavior: 



1. The interaction/bath dynamics is such that there 
is very little relaxation. 
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2. The system as a whole relaxes to some state 
(which may or may not be close to the ground 
state) and this state is a complicated superposi- 
tion of the states of the central system and the 
environment. 

3. The system as a whole relaxes to a state that is 
(to good approximation) a direct product of the 
states of the central system and a superposition of 
states of the environment. In this case there are 
two possibilities: 

(a) The central system does not relax to its 
ground state; 

(b) The central system relaxes to its maximally 
entangled ground state. 

Only case 3b is special: The environment and central 
system are not entangled (to a good approximation) 
but nevertheless the decoherence induces a very strong 
entanglement in the central system. In this paper, we 
demonstrate that, under suitable conditions, dissipation 
free decoherence forces the central system to relax to a 
maximally entangled state which itself, shows very little 
entanglement with the state of the environment. 

Most theoretical investigations of decoherence have 
been carried out for oscillator models of the environ- 
ment for which powerful path-integral techniques can 
be used to treat the environment analytically [2] . On 
the other hand, it has been pointed out that a magnètic 
environment, described by quantum spins, is esscntially 
diffcrent from the oscillator model in many aspeets [7j. 
For the simplest model of a single spin in an external 
magnètic ficld, some analytical results are known J7J- 
For the genèric case of two and more spins, numerical 
simulation |SJ |Ü] is the main source of theoretical infor- 
mation. Not much is known now about which physi- 
cal properties of the environment are important for the 
cfficient selection of pointer states. Recent numerical 
simulations 9 confirm the hypothesis ^ïï| on the rele- 
vance of the chaoticity of the environment but its effect 
is actually not dràstic. 

In this paper, we report on the results of numerical 
simulations of quantum spin systems, demonstrating the 
crucial role of frustrations in the environment on deco- 
herence. In particular, we show that, under appropriate 
conditions, decoherence can cause an initially classical 
state of the central system to evolve into the most ex- 
treme, maximally entangled state. We emphasize that 
we only consider systems in which the total energy is 
conserved such that the decoherence is not due to dissi- 
pation. 



We study a model in which two antiferromagneti- 
cally coupled spins, called the central system, interact 
with an environment of spins. The model is defined by 

H = H c + H e + Hi n t , H c = — JSi ■ S 2 

2V-1 N 
2 N 

^ = -EEE A S )s ^> (i) 

i— 1 j — 1 OL 

where the exchange integrals J < and 0,*$ de- 
termine the strength of the interaction between spins 
S„ = (S%,S%,S*) in the central system (H c ), and the 
spins I„ = (1^,1%, I*) in the environment (H e ), respec- 
tively. The exchange integrals A^ control the interac- 
tion (Hi n t) of the central system with its environment. 
In Eq. I|T)I. the sum over a runs over the x, y and z 
components of spin 1/2 operators. The number of spins 
in the environment is N. 

Initially, the central system is in the spin-up - spin- 
down state and the environment is in its ground state. 
Thus, we write the initial state as |\&(t = 0)) = | Tj 
)|3>o)- The time evolution of the system is obtained by 
solving the time-dependent Schròdinger equation for the 
many-body wave function |í(í)), describing the central 
system plus the environment. The numerical method 
that we use is described in Ref. It conserves the 

energy of the whole system to machine precision. 

By changing the parameters of model JIJ, we explore 
the conditions under which the central system clearly 
shows an evolution from the initial spin-up - spin-down 
state towards the maximally entangled singlet state. We 
consider systems that range from the rotationally in- 
variant Heisenberg case to the extreme case in which 
H e and H int reduce to the Ising model, topologies for 
which the central system couples to two and to all spins 
of the environment, and vàlues of parameters that are 
fixed or are allowed to fluctuate randomly. Illustrative 
results of these calculations are shown in Figs. 1-4. In 
Table I, we present the corresponding numerical data of 
the energy (*(0)|iï|*(0)) = (*(f)|.ff|*(t))) and of the 
two-spin correlation (S x (í) • S 2 (í)) = (*(í)|Si • S 2 |*(í)). 
For comparison, Table I also contains the results of the 
energy E and of the two-spin correlation (Si • 82)0- in 
the ground state of the whole system, as obtained by 
numerical diagonalization of the Hamiltonian Eq.JIJ. 

We monitor the effects of decoherence by computing 
the expectation value (í'(í)ISi • S 2 |\&(í)). The central 
system is in the singlet state if (Si(í) • S 2 (t)) = —3/4, 
that is if (Si(í) • S 2 (í)) reaches its minimum value. We 
also study the time evolution of the concurrence C(í), 
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Table I. Minimum value of the correlation of the central spins and the energy of the whole system (which is conserved), 
as observed during the time evolution corresponding to the curves listed in the first column. The correlations (Si • 82)0 
and the ground state energy Eo of the whole system are obtained by numerical diagonalization of the Hamiltonian 
Eq.Q. 
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-2.365 


-0.738 


-0.735 




Fig. 1 (color online) Left: Time evolution of the correlation (í'(í)ISi • S2 1 í'(í)) of the two spins in the central system. 
Dashed horizontal line at -1/4: Correlation in the initial state ({^(í = 0) |Si • S2|í'(t = 0)) = —1/4); Horizontal line 
at -3/4: Expectation value in the singlet state; (a) Environment containing N = 14 quantum spins; (b) N = 16; (c) 
N = 18. The parameters QÍ™ and AÍ" are uniform random numbers in the range [— 0.15| J|, 0.15| J|]. Right: Time 
evolution of the concurrence C(í) for three different random realizations of a spin glass environment. The parameters are 
uniform random numbers in the range — 0.15|J| < í^"', AÍ°^ < 0.15|J| and the environment contains N = 14 quantum 
spins. The transition from an unentangled state (C(í) = 0) to a nearly fully entangled state (C(í) = 1) is clearly seen. 



which is a convenient measure for the entanglement of 
the spins in the central system ^2j- The concurrence is 
equal to one if the central system is in the singlet state 
and is zero for an unentangled pure state such as the 
spin- up - spin-down state |12| . 



A very extensive search through parameter space 
leads to the following conclusions: 

• The maximum amount of entanglement strongly 
depends on the vàlues of the model parameters 
fií 1 ? and A^. For the case in which there is 



4 



S. Yuan, M. I. Katsnelson, and H. De Raedt 




600 



Fig. 2 (color online) Time evolution of the concurrence 
C(t) for the case of a frustrated antiferromagnetic en- 
vironment. The interactions of the central system and 
the environment are uniform random numbers in the 



range -0.15|J] < A^ 3 < -0.05|J|. The environment 
contains 14 quantum spins, arranged on a triangular lat- 
tice and interacting with nearest neighbors only. The 
nonzero exchange integrals are uniform random num- 
bers in the range -0.55|J| < < — Q.45|«7|. The 
transition from an unentangled state (C(í) = 0) to a 
nearly fully entangled state (C(í) = 1) is evident, as is 
the onset of recurrent behavior due to the finite size of 
the environment. 



(fi 
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Fig. 3 (color online) Time evolution of the correlation 
(Ç(í)|Si ■ S2 1 NE' (í) ) of the two spins in the central sys- 
tem. Environment containing N — 16 quantum spins. 
Dashed horizontal line at -1/4: Correlation in the ini- 
tial state ((*(í = 0)|Si ■ S 2 |*(í = 0)) = -1/4); Hor- 
izontal line at -3/4: Expectation value in the singlet 
state. For all curves (a-f) A|* 3 = A^ = 0, that is 
H int is Ising like. The vàlues of A' z ^ are: (a) random 
-0.0375 \J\ or 0.0375 |J|, (b-e) random -0.075 |J| or 
0.075 |J|, (f) random -0.15 \J\ or 0.15 |J|. The vàlues 
of 0'°j are uniform random numbers in the range: (b) 
[-0.03751 J\, 0.0375| J|], (a,c) [—0.15] J\, 0.15[ J|], (d,f) 
[-0.3|J[,0.3|J|] and (e) [-|J|, |J|]. 



strong decoherence, increasing the size of the envi- 
ronment will enhance the decoherence in the cen- 
tral system (compare the curves of Fig. l(a,b,c) 
and Fig. 4(d,e)). Keeping the size of the environ- 
ment fixed, different realizations of the random 
parameters do not significantly change the results 
for the correlation and concurrence (right panel 
of Fig. 1). However, the range of random vàlues 
ÍV-j and AÏj°j for which maximal entanglement 
can be achieved is narrow, as illustrated in Figs. 3 
and 4. In Fig. 3 we compare results for the same 
type of Hi n t (Ising like) and the same type of H e 
(anisotropic Heisenberg like), but with different 
vàlues of the model parameters. In Fig. 4, we 
present results for different types of Hi nt and H e . 
but for parameters within the same range. 

• Environments that exhibit some form of frustra- 
tion, such as spin glasses or frustrated antiferro- 
magnets, may be very effective in producing a high 
degree of entanglement between the two central 
spins, see Figs. 1-4. 

• Decoherence is most effective if the exchange cou- 
plings between the system and the environment 



are random (in a limited range) and anisotropic, 
see Figs. 3 and 4. 

• The details of the infernal dynamics of the envi- 
ronment affects the maximum amount of entangle- 
ment that can be achieved [ü] , and also affects the 
speed of the initial relaxation (compare the curves 
of Fig. 3(b,c,d,e), Fig. 4(a,d) and Fig. 4(b,c)). 

• For the case in which there is strong decoherence, 
for the same H e and the same type of Hi„ t , de- 
creasing the strengh of Hi nt will reduce the relax- 
ation to the finial state, and the final state comes 
closer to the singlet state (compare the curves of 
Fig. 3(a,c) and Fig. 3(d,f)). 

Earlier simulations for the Ising model in a trans- 
verse field have shown that time-averaged distributions 
of the energies of the central system and environment 
agree with those of the canonical ensemble at some effec- 
tive temperature ^3 El ■ Our results do not contradict 
these findings but show that there are cases in which 
the central system relaxes from a high energy state to its 
ground state while the environment starts in the ground 
state and ends up in state with slightly higher energy. 
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Fig. 4 (color online) Effect of the symmetry of the ex- 
change interactions and A\j on the time evolu- 

tion of the correlation (*(t)|Si • S 2 |*(í)} of the two 
spins in the central system. Dashed horizontal line at 
-1/4: Correlation in the initial state ((*(í = 0) |Si ■ 
S 2 |*(í = 0)) = -1/4); Horizontal line at -3/4: Cor- 
relation in the singlet state; Other lines from top to 
bottom (at t\J\ = 6000): (a) Ising Hi„t with Ising 
H e , N — 14; (b) Heisenberg-like Hint with Ising H e , 
N = 14; (c) Heisenberg-like Hi nt with Heisenberg-like 
H e , N — 14; (d) Ising H int with Heisenberg-like H e , 
N = 14; (e) Same as (d) except that N = 18. We 
use the term Heisenberg-like Hi n t (H e ) to indicate that 
are uniform random numbers in the range 
[— 0.15 \J\ , 0.15 \J\]. Likewise, Ising Hi nt (H e ) means 
that A% y) = (n^' s) = 0), and A$ (íï$) are ran- 
dom -0.075 \J\ or 0.075 |J|. 



As shown in Fig.4(e), this state is extremely robust and 
shows very little fluctuations. 

For the models under consideration, the efïiciency 
of decoherence decreases drastically in the following or- 
der: Spin glass (random long-range interactions of both 
signs); Frustrated antiferromagnet (triangular lattice 
with the nearest-neighbour interactions); Bipartite anti- 
ferromagnet (square lattice with the nearest-neighbour 
interactions); One-dimensional ring with the nearest- 
neighbour antiferromagnetic interactions. This can be 
understood as follows. A change of the state of the cen- 
tral system affects a group of spins in the environmcnt. 
The suppression of off-diagonal elements of the reduced 
density matrix can be much more effective if the group 
of disturbed spins is larger. The state of the central sys- 
tem is the most flexible in the case of a coupling to a spin 
glass for which, in the thermodynamic limit, an infmite 
numbcr of infinitely closed quasi-equilibrium configura- 
tions exist |151 116| . As a result, a very small pertur- 
bation leads to the change of the system os a whole. 
This may be considered as a quantum analog of the 



phcnomcnon of "structural relaxation" in glasses. This 
suggests that frustrated spin systems that are close to 
the glassy state should provide extremely efncient deco- 
herence. 

To conclude, we have demonstrated that frustrations 
and, especially, glassiness of the spin environment result 
in a very strong enhancement of its decohering action on 
the central spin system. Our results convincingly show 
that this enhancement can be so strong that solely due 
to decoherence, a fully disentangled state may evolve 
into a fully entangled state, even if the environment con- 
tains a relatively small numbers of spins. 



1. D. Giulini, E. Joos, C. Kiefer, J. Kupsch, I.-O. Sta- 
matescu, and H.D. Zeh, Decoherence and the Appear- 
ance of a Classical World in Quantum Theory (Springer, 
Berlin, 1996). 

2. W.H. Zurek, Rev. Mod. Phys. 75, 715 (2003). 

3. J. von Neumann, Mathematical Foundations of Quan- 
tum Mechanics (Princeton University Press, Princeton) 
1955. 

4. R.P. Feynman and F.L Vernon, Ann. Phys. (N.Y.), 24, 
118 (1963). 

5. A.J. Leggett, S. Chakravarty, A.T. Dorsey, et al., Rev. 
Mod. Phys., 59, 1 (1987). 

6. J.-P Paz and W.H. Zurek, Phys. Rev. Lett., 82, 5181 
(1999). 

7. N.V. Prokof'ev and P.C.E. Stamp, Rep. Prog. Phys., 
63, 669 (2000) 

8. V.V. Dobrovitski, HA. De Raedt, MI. Katsnelson, et 
al., Phys. Rev. Lett., 90, 210401 (2003). 

9. J. Lages, V.V. Dobrovitski, MI. Katsnelson, et al., 
Phys. Rev. E, 72, 026225 (2005). 

10. W.H. Zurek, Nature 412, 712 (2001). 

11. V.V. Dobrovitski and H.A. De Raedt, Phys. Rev. E, 67, 
056702 (2003). 

12. W.K. Wootters, Phys. Rev. Lett., 80, 2245 (1998). 

13. R.V. Jensen and R. Shankar, Phys. Rev. Lett., 54, 1879 
(1985). 

14. K. Saito, S. Takasue, and S. Miyashita, Phys. Rev. Lett., 
65, 1243 (1996). 

15. K. Binder and A.P. Young, Rev. Mod. Phys., 58 
801(1986). 

16. M. Mezard, G. Parisi, and M.A. Virasoro, Spin Glass 
Theory and Beyond (World Scientihc, Singapore) 1987. 



